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Abstract. The F8Ia supergiant V810 Centauri is part of 
a long-term high-precision photometric monitoring pro- 
gram on long period variables started twenty years ago. 
Time series analysis of this unique set of 500 data points, 
spanning almost fifteen years in the homogeneous Geneva 
photometric system, is presented. Cluster membership, 
physical parameters and evolutionary status of the star 
are reinvestigated. Radial velocity data do not support 
the cluster membership to Stock 14. Ultraviolet and op- 
tical spectrophotometry is combined with optical and in- 
frared photometry to evaluate the physical parameters of 
the yellow supergiant (T e g = 5970 K, Afboi = -8.5, R = 
420 Rq) and of its BO III companion. From theoretical 
stellar evolutionary tracks, an initial mass of ~25 M is 
estimated for V810 Cen, which is actually at the end of 
its first redward evolution. 

V810 Cen is a multi-periodic small amplitude variable 
star, whose amplitudes are variable with time. The pe- 
riod of the main mode, ~156 d, is in agreement with the 
Period-Luminosity-Colour relation for supergiants. This 
mode is most probably the fundamental radial one. Ac- 
cording to the theoretical pulsation periods for the radial 
modes, calculated from a linear non-adiabatic analysis, the 
period of the observed second mode, ~107 d, is much too 
long to correspond to the first radial overtone . Thus, this 
second mode could be a non-radial p-mode. Other tran- 
sient periods are observed, in particular at ~187 d. The 
length of this period suggests a non-radial g-mode. Then, 
the complex variability of V810 Cen could be due to a 
mixing of unstable radial and non-radial p- and g-modes. 
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* Based on observations collected at the Swiss 40 cm and 
70 cm and at the Danish 1.54 m telescopes, at the European 
Southern Observatory (La Silla, Chile) 

** The photometric data are only available in elec- 
tronic form at the CDS via anonym ous ftp to c d- 
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1. Introduction 

V810 Centauri (=HD 101947 =HR 4511, hereafter 
V810 Cen) is a yellow supergiant at low galactic latitude 
(1=295.18, b=-0.64) in the direction of Carina spiral arm. 

The light variability was discovered by Fernic ( 1976 ) 
as a result of cepheid-like s uperg iants photometric sur- 
vey. Eichendorf & Reipurth ( [1979 ) found a low amplitude 
variation (0.12 mag in V) with a "period" of 125 days 
(but only two successive minima were observed). Dean 
(198(3) confirmed the variability of V810 Cen, but not the 
125 d period. In a preliminary analysis of Geneva long- 
term photometric monitoring, Burki (1994) obtained two 
radial modes at 153 and 104 days with low amplitudes 
(0.07 and 0.04 mag respectively). However, the high resid- 
uals suggested that the amplitudes were variable and/or 
that additional frequencies may be present. 

The pulsation interpretation of the variability was con- 
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troversial. Indeed, Bidelman et al. (1963) suspected the 
star to be a spectroscopic binary and this was confirmed 
by van Genderen ( |1980[ ) who proposed a B spectral type 
companion to account for the UV excess of the GO la star. 
Therefore, the light variation may come from the com- 
panion variability rather than the GO star. This issue was 
settled with IUE spectra from Parsons ( |1981| ). From con- 
tinuum flux fitting and photometric data he found that the 
G supergiant should be 3.2 mag brighter (in the V band) 
than its B companion, thus comforting the cepheid-like 
pulsation hypothesis for V810 Cen. However, the spectral 
type of the B companion was still ambiguous: C IV and 
Si IV absorption lines lead to a B0-B1 lab-lb star while 
continuum flux fitting requires a B0.5-B1 III star (see Sec- 
tion 3). Since the G star luminosity relies on its companion 
bolomctric magnitude and their magnitude difference, the 
variable star luminosity is also ambiguous. 

The perspective of using V810 Cen as a prime candi- 
date to understand the variability of the massive stars was 
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strengthened by the open cluster membership. Moffat & 
Vogt ( [1975D suggested that V810 Cen is a member of the 
open cluster Stock 14 for which Peterson & FitzGerald 
( |1988D derived < E B -v >= 0.26 and V - M v = 12.14. 
But the intrinsic luminosity, derived from the cluster dis- 
tance, is 0.5 mag fainter when compared to the spectro- 
scopic estimates from IUE data. 

In the present paper we re-investigate various aspects 
of this long-period variable star. Section 2 is devoted to the 
discussion of the membership to the open cluster Stock 14 
based on the cluster members radial velocity. IUE final 
archive and visual spectrophotometry are used in Section 4 
to constrain the luminosity and effective temperature of 
the GO supergiant. The variability of V810 Cen is dis- 
cussed thoroughly on the basis of the Geneva photometric 
data, collected over 15 years (Sections 5 to 8). Finally, a 
discussion of the pulsation modes is given in Section 9. 

2. V810 Cen and Stock 14 



the velocity could be variable) for HD 101838 (Feast & 
Thackeray, 1963). A mean velocity of — 6 ± 2 km/s can be 
adopted for the cluster in agreement with the estimation of 



Moffat & Vogt (P-975D postulated that V810 Cen is a mem- 
ber of the open cluster Stock 14 (C 1141-622). This would 
strongly constrain the physical parameters of the super- 
giant. Unfortunately, as we shall see hereafter, this is prob- 
ably not the case. 

With the UBV photometry from Moffa t fc Vo gt ( |l975| ), 
Turner ( |l982j) a nd Pe terson & FitzGerald ( |l98Sj ), collected 
by Mermilliod ( |l997 ) in his Cluster Database, an d usin g 
the stellar evolutionary tracks from Schaller et al. ( 1992| ), 
the following parameters can be derived for Stock 14 : 
E{B -V) = 0.26, distance= 2.63 kpc, age= 1.4 • 10 7 yr. 
These values are consistent with those found by Moffat 
& Vogt ( p975| ), FitzGerald k Mill er (|l98|), Lynga (^98^) 
and Peterson & FitzGerald ( [19881 ). 

Assuming cluster membership, the absolute magnitude 
of V810 Cen (both components together) would be My = 
-7.88, with an UBV intrinsic color of (B - V) = 0.526. 



According to the analysis of van Genderen (1981), the blue 
component has an absolute magnitude of My = —4.25 and 
an intrinsic color index of (B — V)o = —0.28 (09-BO type 
star), thus the red supergiant would have My = —7.84 
and (B — V)o — 0.56, corresponding to a GOIa type star. 
With these values, the blue component would be close to 
the cluster turnoff of the main sequence, whereas the red 
supergiant remains close to the original position of the 
double system, i.e. close to the blue loop of the evolution- 
ary track, in the core helium burning phase. 

As we shall see from the analysis of the radial veloc- 
ity data, it is probable that V810 Cen is not a member of 
Stock 14. The radial velocities of 4 bright B-type members 
of this cluster are known: -8 km/s (1 measurement with 
an uncertainty of at least 5 km/s) for HD 101994 (Bus- 
combe & Kennedy, 1969| ), — 6.4±1.9km/s for the eclipsing 
SB2 s ystem V346 Cen, HD 101897 (Hernandez & Sahade 
( |1978|) , -3 ± 12 k m/s (4 measurements) for HD 101964 
(Feast et al., 1957) and —10 ± 21 km/s (5 measurements, 



Lynga (1987). This value relies strongly on the 7 velocity 
of the SB2 system V346 Cen from 44 V r measurements. 

V810 Cen has 90 radial velocity data points between 
HJD 2 444621 and 2 449 915 from Coravel spectropho- 
tometer attached to the 1.54 m Danish telescope in La 
Silla. The uncertainty on each measurement is about 0.4 
km/s, whereas the dispersion of the data is much larger, 
due to the complex pulsation of V810 Cen, which is the 
only component of the system measured with Coravel. 
The mean velocity is —16.7 ± 2.6 (s.d.). For further cal- 
culation, the annual means of the velocity will be used. 
The values of 13 annual means range from -13.2 km/s to 
-20.8 km/s, with standard deviations between 0.5 and 3.5 
km/s. 

Older measurements of the radial velocity of V810 Cen 
are: +11.7 km/s in 1908 (2 measurements) and +6. 9 km/ s 
in 1911 (1 measurement) from Campbell & Moore ( 1928 ), 
and -17.4 km/s in 1946 (10 measurements), -14.4 km/s (3 
measurements) in 1947 and -1 0.8 km /s (2 measurements) 
in 1959 from Bidclman et al. ( 1963 ). 

If all the radial velocity data are taken into account, 
the possible orbits with a 7 velocity corresponding 
to the velocity of the cluster members (7 = —6 
km/s, see Fig. Q) give values for the mass function 
/(to) = (02 sin(i)) 3 /P 2 = (toi sin(z)) 3 /(mi +to 2 ) 2 which 
are far too large: the values for mi (B-type component) 
are larger than 60 Mq, if we assume an to 2 of about 25 
M Q . 

If the data obtained in 1908 and 1911 are excluded (these 
are 20.3 A /mm spectra taken on photographic plates), 
the remaining velocities are compatible with a constant 
velocity (taking into account the pulsation), or with 
a small amplitude orbit. But the 7 velocity would be 
close to -16 km/s instead of -6 km/s (cluster velocity). 
Furthermore, it is noteworthy that our recent data have 
a mean value (—16.7 ± 2.6 km/s) in good agreement 
with the mean value of the data obtained in 1946 and 
1947 (—15.8 ± 4.2 km/s), suggesting that the systematic 
velocity of the system is indeed 10 km/s lower than the 
mean velocity of the cluster. 

In conclusion, our radial velocity analysis allows to ex- 
clude the membership of V810 Cen to the cluster Stock 14. 



3. The blue companion 

The existence of a blue companion was suggested by Par- 
sons & Peytremann ( 1973 ) and van Genderen ( 1980 ) from 
photometric data and conclusively proved by Eichendorf 
et al. ( 1981 ) and Parsons ( 1981 ) on the basis of IUE spec- 
tra. 

From Si IV and C IV line profiles and intensity anal- 
ysis, these authors conclude that the spectral type of this 
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Fig. 1. Radial velocities from 1908 to 1995 collected from 
various sources. Only mean values are shown. The impor- 
tant scatter in the last set of data is due to V810 Cen 
pulsation. The fitted orbit (full curve) would agree with 
the 1908-1928 data but is physically unlikely (see text). 



blue companion is B0-B1 lab-lb, i.e. a star of absolute 
magnitude My ~ —6. Furthermore, Parsons (1981) fitted 



synthetic spectrum on IUE data and UBVRIJKL photom- 
etry and determined a V magnitude difference between the 
two components of AV ~ 3.2 (the blue component being 
fainter). As noted by Parsons, there is an inconsistency 
between V810 Cen membership to the cluster Stock 14 
and the luminosity class of the companion. Indeed, the 
absolute magnitude of the two components would be -7.8 
(V810 Cen) and -4.6 (blue companion) if the membership 
is accepted and the value AV = 3.2 taken into account. 
In that case, the luminosity of the blue companion corre- 
sponds to a giant and not to a supergiant (as it is indicated 
by Si IV and C IV lines). I n orde r to reconcile spectral and 
cluster luminosity, Turner ( 1982 ) suggested that the B star 
stellar wind could be abnormally strong for its luminosity, 
as it is the case for t Sco (BO V) and £ Oph (09 V). If the 
same phenomenon applies here, this star could be a giant, 
with My ~ —4.6, in agreement with the membership to 
Stock 14. 

However, as noted in the previous section, the ra- 
dial velocity analysis do not support the membership to 
Stock 14. Thus, it appears to us that the most logical solu- 
tion is to accept : i) the luminosity class of the blue com- 
panion, probably a giant BO star (see section 4), and its 
absolute magnitude My ~ —5.1 (Schmidt-Kaler, 1981); 



Fig. 2. Reddened and scaled Kurucz models (dashed lines) 
assuming a BO III hot companion. The full line is the sum 
of the BO III and F8Ia Kurucz model. IUE and Kiehling's 
data (crosses), the flux derived from Geneva photometry 
(filled squares) and IRAS flux (open squares in the insert) 
are also given. The error bars are the upper and lower flux 
observed in Geneva photometry whereas for IRAS data 
they represent the estimated error. The arrows indicate 
upper limits for the flux. 



the absolute magnitude of V810 Cen is My ~ 
the star is located behind the cluster Stock 14. 



-8.4 and 



4. Physical parameters of the two components 

Theoretical energy distributions for each of the two com- 
ponents, t aken f rom the grid of stellar atmosphere models 
of Kurucz ( 1994 ), have been fitted to the observed UV, op- 
tical and IR spectrophotometric data of V810 Cen (both 
components together). 

In the UV domain, three low dispersion, large aper- 
ture, spectra obtained with IUE satellite and reduced with 
NEWSIPS (Nichols & Linsky, |l99l| ) are available in the 
IUE Final Archive (obtained from NASA Data Archive 
and Distribution Service), namely lwp6460, lwp28437 and 
swp26455. The two lwp spectra have been averaged in the 
2122-3348 A range and the strong chromospheric Si IV 
and C IV lines in the swp spectrum have not been used 
since they cannot be well described by the standard Ku- 
rucz model atmosphere. 

In the optical domain, a 10 A resolution spectrum is 
available in the range 3250 to 8600 A (Kiehling, |1987|). 



ii) the V magnitude difference between the two compo- 
nents, AV ~ 3.3 (see next section). With these values, 



These data are accessible through the Strasbourg Stellar 
Data Center. Points flaged as uncertain by the author have 
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Star 


Sp. type 


E(B - V) 


T c {{ 


T c ff 


AlogTeff 








Evans & Teays 


This paper 




a Per 


F5Iab 


0.04 


6 270 


6 430 


0.011 


(3 Aqr 


GOIb 


0.03 


5 560 


5 660 


0.008 


9 Peg 


G5Ib 


0.13 


5110 


5170 


0.005 



Table 1. Temperature estimates for 3 standard supergiants compared to the values derived by Evans & Teays (1996). 
Our temperatures are, in average, 100 K higher than Evans & Teays estimates. 



Sp T efr ,i f logsi f Ci r effi2 log 52 C 2 x 2 



09 V 


33 000 


4.0 


20.286 


6 070 


1.2 


17.204 


2983 


B0 V 


30 000 


4.0 


20.160 


6 040 


1.0 


17.198 


3087 


Bl V 


25 400 


4.0 


19.893 


5 960 


0.4 


17.182 


3654 


09 III 


32 000 


4.0 


20.249 


6 060 


1.1 


17.202 


2975 


B0 III 


29 000 


3.5 


20.119 


6 010 


0.7 


17.190 


2539 


Bl III 


24 000 


3.5 


19.812 


5 905 


0.0 


17.171 


4236 


09 lab 


32 600 


4.0 


20.271 


6 060 


1.1 


17.203 


2979 


B0 lab 


25 000 


3.0 


19.898 


5 910 


0.0 


17.171 


4015 


Bl lab 


20 800 


3.0 


19.587 


5 810 


0.0 


17.151 


9278 


B2 lab 


18 500 


3.0 


19.387 


5 760 


0.0 


17.143 


12405 



Table 2. Parameters for various fits to IUE and visible spectra. The microturbulence is set to 4 km/s and solar 
metallicity is assumed. The color excess is fixed to E(B — V) = 0.26. Fixed parameters are quoted with an f. The best 
fit is achieved for a B0 III and F8 la couple. 




Fig. 3. V810 Cen position in the HR diagram assuming 
a B0 III hot companion. The "error bars" corresponds to 
250 °K and 0.5 mag. Evolutionary tracks are from Meynet 
et al. (1994) and the instability strip (dashed lines) is from 
Chiosi et al. (1993). Linear non-adiabatic stability analysis 
has been performed on nine models (diamonds) on the 25 
Mq track (see also Table 0). Dotted lines are the iso- 
periods in the 50 to 200 days range (25 days step) for the 
fundamental mode calculated from Schallcr (|199C| ). 



not been used. In the IR domain, the photometric data 
from IRAS satellite (IRAS Points Source Catalog, |l988| ) 
are reliable only at 12 fim, whereas the fluxes at 25, 60 and 
100 fim are only upper limits. In addition, the photometric 
measurements in the 7 colors Geneva system have been 
used to estimate the effect of the variability. 



No detailed study of the interstellar extinction in the 
UV domain is available towards V810 Cen. Assuming that 
the star is not too far away from the open cluster Stock 14, 
which is at a di stanc e of 2.7±0.2 kpc according to Peterson 
& FitzGcrald ( 1988), the mean Galactic extinction law of 
Krelowski & Papaj (|l992|) was adopted. This law is suited 
for stars at 2 to 6 kpc from the sun. 



For each absolutely calibrated flux of the extracted 
IUE spectra an internal error estim ate is given in the 
MXLO files (Nichols et al., |l993|) . The spectra have 



been co rrecte d for the systematic deviations described by 
Bohlin ( 1996 ). From his Fig. 1, the following values for the 
IUE external error have been adopted: 3 %, 5 %, 2 %, and 
7 % respectively in the ranges 1100-1900 A, 1900-2400 A, 
2400-3100 A and 3100-3300 A . For the optical spectrum, 
Kiehling's external and internal error estimates have been 
adopted. 

To reproduce the observed composite spectrum, with 
two model atmospheres and the mean interstellar extinc- 
tion law, seven parameters should be estimated: the effec- 
tive temperature, gravity and angular diameter for each 
star and the color excess. The estimate has been per- 
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formed through a standard least squares procedure mini- 
mizing the following quantity: 



[log(FP) - log(Ff )] 5 



(1) 



where and Ff 1 arc the unreddened observed and theo- 
retical fluxes at the wavelength Ai, and <7i is the error esti- 
mates on log(F° fcs ), where F° bs is the observed, reddened, 
flux. To compare the observed and theoretical fluxes the 
former has been resampled at the model frequencies. The 
theoretical flux is given by: 



F. 



Hi 



F™° d {T ugi ) 10- 



F™° d {T 2 ,g 2 ) 10 



-c 2 



where F™ {T,g) is the flux from Kurucz's atmosphere 
models at A; for given values of T g and logg. The indices 
1 and 2 refer respectively to the blue and red (V810 Cen) 
components. The Kurucz's models used are those with the 
solar abundance and a microturbulence of 4 km/s, well 
suited for supergiants. The parameter C is related to the 
stellar angular diameter 9 through: 



9 = 2- 2.06 10 8 10 



-0.5 C 



milliarcsec 



The unreddened observed flux is given by: 
log(F?) = log(^ ofcs ) + 0.4 At 



(2) 



(3) 



where Ai is the interstellar extinction in magnitude, calcu- 
lated from the extinction law of Krelowski & Papaj ( 1992Q , 
assuming R = A V /E(B -V) = 3.1. 

To estimate the accuracy of the present method, 
we used the spectrophotometric data from Glushneva 
et al. fll992| ) for 3 supergiant stars: a Per, (3 Aqr and 
9 Peg. Their parameters have been derived by Evans & 
Teays ( |l996|) on the basis of ultraviolet measurements 
(IUE data), optical and infrared photometry (BVRIJHK 
bands). We have adopted the same values for the gravity 
and the microturbulence as Evans et al. (logg = 1.5 and £ 
=4 km/s for the 3 stars). Table [j] shows that our temper- 
atures are slightly larger than Evans et al. values by 100 K 
on the average. However, the mean difference is small, less 
than 0.01 in logT c ff, i.e. less than one spectral subclass. 
We adopt that value of 100 K for the typical uncertainty 
of our temperature determinations of V810 Cen. 

According to Eichendorf & Reipurth (197£) and 



Turner fll982j ) the blue component of V810 Cen is in the 
ranges B0-B1 in spectral type and V-Iab in luminosity 
class. The parameters T e s,i and loggi of the blue com- 
ponent were fixed, and the minimization of the x 2 value 
was performed with the free parameters Ci, T G ff, 2, log 52 
and C2. The computed ranges in spectral type and lu- 
minosity class for the blue component were 09-B2 and 
V-Iab. In a first step, the color excess was fixed to 0.26, 



the mean value for the cluster Stock 14. Table g shows 
that the best fit is achieved for a blue component of type 
B0 III. The corresponding values for the red component 
are T ff,2 = 6 010K and log 52 = 0.7, corresponding to an 
F8 supergiant. 

In a second step, the color excess was varied from 0.20 
to 0.30, for the same parameters of the blue component, 
i.e. those of a B0 III star. As shown by the % 2 values in 
Table ||, the fitting procedure is weakly dependent on the 
color excess value : any value of E(B — V) between about 
0.20 and 0.28 can be accepted. 

From the relation : 



log{R 2 /R x ) = 0.S{C l -C 2 ) 



(4) 



we derive the radius ratio R2/R1 =31.1 which clearly 
indicates that the luminosity class of the red component 
is la. Thus, the best solution for the two components of 
V810 Cen is : 

Blue component : B0 III, T cS = 29 000 ± 1 000 K 
Red component : F8Ia, T cff = 5 970 ± 100 K 

This solution is plotted in Fig. ^ where reddened model 
fluxes (dashed lines) and their sum (full line) are com- 
pared to the fluxes in the 7 Geneva passba nds ( full boxes, 
computed according to Rufener & Nicolet, 1988| ). Further- 
more, IRAS photometry is available for the point source 
11410-6212 associated to V810 Cen and it is compared to 
the model in the insert. The error bars, associated with 
Geneva photometry, are the maximum and minimum ob- 
served fluxes, while those plotted at 12 /im are the 16 % 
uncertainties quoted in the IRAS point source catalogue. 
Those data marked by an arrow are upper limits fluxes. 

The absolute and bolometric magnitudes have been 



derived from the spectral type (Lang, 1992) for the blue 
component and using Ay = 3.3 for the red component 
(V810 Cen). The values of the stellar radii have been 
calculated from T c g and M\, Q \. Thus we have : 



Blue comp. : My — —5.1, Afboi = 
Red comp. : My = —8.4, Mboi = 



-8.0, R= 14 R Q 
-8.5,i? = 420 R 



Note that the ratio R2/R1 is 29.8, very close to the 
value derived above from the values of the parameter C. 
The two components can then be placed in the logL vs. 
log T ff diagram (see Fig. ||) together with the evolutionary 
tracks from Meynet et al. (1994). Initial masses of 25 ± 5 
Mq can be derived for both components. 



As already noted by Eichendorf & Reipurth ( |1979[ ), 
V810 Cen is located near the blue edge of the instability 
strip (in Fig. || the limits of the strip for th e fun damental 
mode, dashed lines, are from Chiosi et al. ( |1993; )). 

With the values of the absolute magnitude of both 
components and of the mean visual apparent magnitude 
(V = 5.021), we derive for V810 Cen a distance of 3.3 
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5. The photometric data 

V810 Cen was measured 512 times in the seven 
ters of t he G eneva Photometric system (Golay. 



fil- 



1980 




Rufener, |1988| ) from JD 2 444 544 (Nov. 1980) to 2 450 051 
(Dec. 1995) with a lack of observations from 2 446 638 to 
2 447861 (see Fig. |J). These data come from the succes- 
sive Swiss telescopes (40 cm and 70 cm) at ESO La Silla 
Observatory (Chile), equipped w ith th e P7 photoelectric 
photometer (Burnet & Rufener, 1979). The data reduc- 



tion has been made according to the method described 
by Rufener ( 1988j ). Only the 499 measurements having 
weights in magnitude (Q) and in colors (P) larger than 
have been used. 

The photometric data are available in electronic 
format at the CDS via anonymous ftp to cd- 



sarcu.strasbg.fr (130.79.128.5) or via bttp://cdsweb.u- 



strasbg.fr/Abstract.html 



6000 8000 
HJD a'440'000 



Fig. 4. Geneva photometric data for V and [B-V]. The 
dashed line in figure a shows the long-term brightening of 
0.0091 mag/y. 



to 3.5 kpc (for E(B-V) in the range 0.28 to 0.24). Thus 
V810 Cen is located 0.6 to 0.8 kpc behind the stellar clus- 
ter Stock 14. Is this location in contradiction with the color 
excess values (0.26 for Stock 14, 0.20-0.28 for V810 Cen) ? 
The answer is no, because Stock 14 is non-uniformly red- 
dened. For the member stars of the cluster, Peterson & 
FitzGerald ( 1988j ) found a standard deviation of 0.021 on 
the individual values of E(B — V). This value corresponds 
to a 3(7 interval in E(B - V") of 0.19-0.32. Moreover, the 
two stars of Stock 14 which have the lowest reddening, 
E(B - 1Q =0.18 (star 26) and 0.21 (star 27) according 
to Turner ( ^98^ ), are located in the sky very close to 
V810 Cen (stars 33, 34 and 35, are also close, but cannot 
be used because of their unreliable photometric data). 



E(B-V) 


x 2 


Ci 


Teff,2 


log (#2) 


c 2 


0.20 


2795 


20.314 


5 890 


1.1 


17.216 


0.22 


2581 


20.249 


5 930 


1.0 


17.207 


0.24 


2495 


20.184 


5 970 


0.9 


17.199 


0.26 


2639 


20.119 


6 010 


0.7 


17.190 


0.28 


2721 


20.054 


6 040 


0.6 


17.181 


0.30 


3069 


19.988 


6 070 


0.5 


17.171 



Table 3. Variation of the \ 2 with respect to E(B — V) if a 
B0 III companion is assumed (see Table ||) . The microtur- 
bulence is set to 4 km/s and solar metallicity is assumed. 



Table |4] displays the mean magnitude values and rms 
for the 7 filters for V810 Cen and 3 standard stars. The 
corrected rms value <7o for V810 Cen is given by <7q = a 2 — 
[(l/3)(cr s i +<t S 2 +ct S 3)] 2 , where er is the rms of V810 Cen 
((To differ from a by less than 0.001). At this stage, the 
essential characteristics of the variability of V810 Cen is : 
i) an increase of ao, thus of the global amplitude from G 
to Bl filters; ii) an amplitude in U smaller than in Bl, 
B and B2. This apparent decrease of the amplitude in U 
is due to the constant flux from the B-type component, 
which contribute significantly to the total flux only in this 
filter. 

Our monitoring of V810 Cen can be characterized as 
follows : 

- Three measurements have been obtained in Febru- 
ary and March 1976, and in February 1977 (HJD from 
2 442 817 to 2 443182). 

- From October 1980 to July 1986 (HJD from 2 444 544 
to 2446 639), the air mass of the measurements was re- 
stricted to values smaller than about 1.7. Thus the star 
was measured each year, only between November and July. 
On the average, 25 measurements per year were obtained. 

- From December 1989 to October 1991 (HJD from 
2 447860 to 2 448 559), a peculiar observational effort was 
done in order to have a monitoring as continuous as pos- 
sible (the star is observable during the whole year from 
La Silla). On the average, 9 measurements per month were 
obtained. These data are shown in Fig. ^a. Variations with 
a characteristic time in the 100-150 days range are clearly 
exhibited. In addition, the amplitude is variable. 

- From November 1991 to December 1996 (HJD from 
2 448 590 to 2 450 051), the monitoring was less intensive. 
Due to several periods without data, the average number 
of measurements per month is 3. 
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Table 4. Weighted mean value and rms of the 499 measurements of the 7 magnitudes of V810 Cen in the Geneva 
Photometric System. The rms is also given for 3 standard stars measured during the same nights than V810 Cen. 





Filter 


u 


Bl 


B 


B2 


VI 


V 


G 




Ao 


[A] 


3464 


4015 


4227 


4476 


5395 


5488 


5807 


V810 Cen 


Mean mag. 


6.797 


6.155 


5.035 


6.323 


5.785 


5.021 


6.005 


V810 Cen 


corrected rms (ao) 


0.055 


0.081 


0.070 


0.062 


0.046 


0.044 


0.040 


HD 93540 


rms 


(<T s i) 


0.0070 


0.0061 


0.0057 


0.0058 


0.0050 


0.0043 


0.0048 


HD 94510 


rms 


(<T S 2) 


0.0083 


0.0065 


0.0060 


0.0058 


0.0048 


0.0040 


0.0048 


HD 102350 


rms 


Os3) 


0.0068 


0.0056 


0.0047 


0.0049 


0.0044 


0.0035 


0.0046 



6. The photometric variability 



remained constant. 



Three methods have been used to analyze the variability 
of V810 Cen : 

- The Date Compensated Discret e Fou rier Transform 
DCDFT described by Ferraz-Mello (|l98l|) . This method 
avoids the missmeasurement of the amplitudes and 
mean magnitude encountered whe n the classical Discrete 
Fourier Transform (Dee ming, |l975) is applied to unevenly 
spaced data (see Foster, 1995| , 1996a ). 



- The CLEANEST method of Foster ( |l995| ). The 

spurious peaks in the Fourier analysis induced by the data 
sampling are eliminated by an iterative deconvolution 
process. The iteration has been stopped when the residual 
peaks in DCDFT were below the 1% confidence level 



limit according to Foster (1996a) 



The Weighted Wavelet Z-transform method 



(WWZ) of Foster ( |l996b| ). Schematically, the data are 
weightened with a gaussian function (with a ~ 200 
d to 300 d in our present case) centered at a given 
time t. The Fourier Transform is then performed for 
successive t values and the WWZ maximum (WWZ max ) 
indicates the dominating frequency ^(WWZ max ) (noted 
hereafter ^ max )- The amplitude of ^ max is estimated by 
the means of the Weighted Wavelet Amplitude statistics 
(WWA max ). The WWZ method allows to analyze the 
changes in amplitude or/ and frequency of the light curve. 

In order to check the period and amplitude variabil- 
ities, the data have been splited into two samples: sets 

1 and 2 span the ranges in HJD 2 444544-2446 639 and 

2 447860-2450 051 respectively (see Fig. |). A careful ex- 
amination of the data shows that the observed variations 
cannot be simply described as a stable, multiperiodic 
function. The main characteristics of the variability of 
V810 Cen during our survey are the following : 

- Long-term variation : the mean luminosity of the 
star increased regularly from February 1976 to July 1986 
(see Fig. ^). The rate of brightening was roughly 0.009 
mag/year. Since November 1991, the mean luminosity 



- Modes at ^156 and ^107 days : Figs. ||a and 
b present the DCDFT for both sets. We note two main 
peaks at frequencies 0.0064 and 0.0093 d _1 . The fre- 
quency difference is close the 1 y^ 1 alias (~ 0.00274 d _1 ), 
suggesting that one of these may be an alias. However 
both of them are real because they are kept by the 
CLEANEST iterative method. 

- Other modes : the CLEANEST spectra (Figs. |c 
and d) reveal a third period at 115 d for set 1, and five 
new periods at 89, 129, 167, 185 and 234 days for set 2. 
These peaks appear thanks to the better time sampling 
of the second set. 

A curve of the form : 



f(t) =Y^ A i COS[2-KVi{t - to) + fa] 



(5) 



has been fitted to the observations in sets 1 and 2 and 
the resulting amplitudes are given in Table |. The 
residual standard deviation is 0.019 mag for both sets, 
a quite large value compared to the accuracy of our 
measurements (0.004 mag, see Table @). As we shall see 
in the following, this is due to amplitude or/and period 
variation. 



7. Analysis of the mode variations 

The data density is large enough in set 2 (see Fig. |]a) 
to allow an analysis of the variations of the frequencies 
and/or amplitudes of the dominant modes. Two kinds of 
analysis have been done : 

- DCDFT method is performed in five successive inter- 
vals. In each of them (see Fig. ^,), the three most impor- 
tant modes have been identified and their frequencies and 
amplitudes are plotted as full symbols in Figs, and c 
(the decreasing amplitude are symbolized with triangles, 
squares and dots respectively). 

- WWZ method described previously is applied to 
the entire interval. In our case, the WWZ method gives 
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Table 5. Frequency, period and amplitude estimate from 
CLEANEST for set 1 (HJD 2 444 544-2446 639) and 2 
(HJD 2 447860-2450 051). 





Frequency [c/d] 


Period [d] 


Amplitude [mag.] 


E 


Set 1 


0.00641±2e-05 


156.0±0.5 


0.028±0.002 






0.00964±2e-05 


103.7±0.2 


0.028±0.002 






0.00873±3e-05 


114.6±0.4 


0.019±0.002 


0.012 


Set 2 


0.00668±2e-05 


149.8±0.2 


0.036±0.002 






0.00912±3e-05 


109.6±0.2 


0.024±0.002 


5. 0.01 




0.01119±3e-05 


89.4±0.2 


0.020±0.002 


• 0.00B - 




0.00541±3e-05 


184.8±0.5 


0.024±0.002 


0.006 - 




0.00597±2e-05 


167.4±0.4 


0.025±0.002 






0.00428±4e-05 


233.5±1.3 


0.015±0.002 






0.00777±5e-05 


128.7±0.5 


0.013±0.002 


— 0.06 



a 

SET 1 



DCDFT 




CLEANEST 



3ET 1 



b 

SET 2 



DCDFT 




fVwy 



CLEANEST 




0.005 0.01 0.015 0.020 0.005 0.01 0.015 0.02 
v[c/d] v[c/d] 



Fig. 5. a and b: Date-compensated discrete Fourier trans- 
form (DCDFT) amplitudes for sets 1 and 2. c and d: 
CLEANEST for sets 1 and 2, it is a combination of the 
subtracted periods and amplitudes (vertical bars each la- 
beled with their period) and the DCDFT of the residu- 
als (full line). CLEANEST has been applied until all the 
peaks of the DCDFT power are below the 1 % confidence 



level limits (computed according to Foster, 1996a 



reliable results only for the mode with the largest power 
(WWZ max ). The variations of the dominating frequency 
i/ max and its amplitude WWA max are shown as solid lines 
in Figs, ^b and c. 

The main points revealed by this analysis are : 

1. The mode at ~107 d (~ 0.0093 d" 1 ) is the most im- 
portant one before HJD 2 448400, while the mode at 
~156 d (~ 0.0064 d _1 ) dominates afterwards (with 



t " "s* ■ ft 





8000 8500 9000 

HJD 2'440'000 



Fig. 6. a: V-band photometry for data set 2, with sub- 
set division (vertical dotted lines), b and c: dominating 
frequencies (WWZ maximum) for a given "mean time" r 
and its corresponding amplitude WWA (full lines). Trian- 
gle, square and dot symbols (sorted according to decreas- 
ing amplitude) are the tree main frequencies found with 
DCFT analysis for each subset shown in a. 

short exceptions around HJD 2 449 500, 2 449 580 and 
2 449 950) . At the end of our survey, it appeared that 
the mode at ^107 d was again the most significant one 
(see Fig. |b). 

2. The amplitude of the most important mode varies 
almost continuously. However, a stable amplitude at 
~0.06 mag was observed in the HJD interval 2 448 700- 
2 449 100, for the mode at ~156 d (see Fig. |c). 

3. The DCDFT analysis (Fig. |b) shows that a mode at 
- 0.0045 - 0.0055 d" 1 (period between 225 and 180 d) 
is always present, but is never the most important one. 

In conclusion, V810 Cen is a multi-periodic small 
amplitude variable star, whose periods and ampli- 
tudes are variable with time. The periods of the two 
dominant modes are ^156 and ~107 d. The 7 years inten- 
sive survey (1989-1996) also reveals a third mode, with a 
period in the range 180-225 d. 

8. The light and color curve parameters 

The amplitude ratios and the differences in phase between 
the light V and color [B — V] curves are interesting for the 
variability description, in particular for the mode identifi- 
cation. Because of the multiperiodic character of the light 
curve and of the amplitude variation described in the pre- 
vious sections, the amplitude and phase of the principal 
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modes have been fitted in various time intervals through 
equation (5), setting the number of frequencies (n) to 3. 

According to Fig. [| the parameters of the mode at 
~107, ~156 and ~187 d have been calculated respectively 
in the time intervals in HJD 2 447 850-2 448 400, 2 448 400- 

2 449 400 and 2 447 850-2 450 100. The results are given in 
Table ^. The light and color curves for the two principal 
modes at ~107 and ~156 days are shown in Fig. ^. As we 
see, the light and color variations are clearly in phase, i.e. 
the star is bluer when brighter. This is an indication of 
cepheid-like pulsation nature. However, it must be noted 
that, taken into account the uncertainties on A<p in Table 
||, the values of the phase difference <fiv — <t>B-v can be 
slightly negative, zero or slightly positive. According to 
Balona & Stobie ( 1979| , 1980Q , this indicates that these 
two modes can be radial (A</> < 0), but could also be non- 
radial, with odd (Atp = 0) or even (A(f> > 0) values of the 
spherical harmonic order 

The third mode, at ^187 days, shows a clearly positive 
value of Acj). This might indicate a non-radial quadrupole 
mode (1 = 2). However, due to the amplitude and period 
changes in the modes of V810 Cen, the mode identification 
must not be made on the basis of Atfi values only (see next 
section) . 

When the time interval is short enough and the mon- 
itoring sufficiently dense, very good fitted light and color 
curves can be achieved. This is illustrated in Fig. |[ for 
the interval in HJD 2 447850-2448150 (the first interval 
in Fig. ||a). The solid curve in Fig. || represents a fit with 

3 modes, the most important one being at ~ 107 d. The 
residual standard deviation of the observed values around 
the fitted curves are 0.0050 mag in V and 0.0046 mag 
in [B — V}. Thus, in that case, the three modes describe 
completely the variability of V810 Cen (see Table 0). 



9. The origin of the variability 

9.1. Comparison with supergiants and cepheids 
calibrations 

— The Period Luminosity Colour PLC relations. 

The largest amplitude period (~156 d) is in very good 
agreement with the empirical PLC relations for yel- 
low supergiants. Taking the parameters derived in Sec- 
tion 4, Mboi = -8.5 an d T cff = 5970 K, the relation 
from Maeder & Rufener fll972| ), logP = -0.346M bo i- 
3 log T c ff + 10.60, gives P = 163 d. Furthermore, as- 
suming an evolutionary mass of ~20 M (see the 



end of this Section), the relation from Burki ( 1978| ), 
logP = -0.38M bo i -3 log T efi -0.5 log M+ 10.93, gives 
P = 152 d. 

The agreement is not so good with the Period- 
Luminosity relation for galactic cepheid stars (Gieren 
et al., |199^| ) : M v = -1.371 - 2.986 log P (scatter : 
0.26). Using the My value obtained in Section 4, 
the predicted period is 226 d. This discrepancy with 




b,- 





T 

so 

E 




0.5 1 



0.5 1 
Phase 



Fig. 7. Light and color curves of the two main modes 
at 107.4 d (time interval: HJD 2 447850-2448400) and 
156.4 d (time interval: HJD 2 448400-2 449 400). 



the observed period (~156 d) is quite normal since 
V810 Cen is not located in the center of the cepheid 
instability strip, but rather in its blue border. For that 
reason, a color term is necessary in the PLC relation 
for supergiants (see Burki, 197§| ). Thus, the observed 
main mode is in agreement with the observa- 
tions concerning pulsating supergiants. 
The Period Radius relation for cepheids. Gieren 
et al. ( 1989] ) have established such a relation on the 
basis of 101 classical cepheids in the period range 3- 
45 d, studied from the visual brightness technique : 
logi? = 1.108 + 0.743 logP (scatter a = 0.071). It is 
noteworthy that, despite the very large value of its 
main period (~156 d), which is well outside the period 
range covered by the classical Galactic cepheids, the 
radius of V810 Cen found in Section 4 (logi? = 2.62) 
is in agreement with this global relation for cepheids 
(logi? = 2.74), the difference being only 1.6 a. This 
agreement is a strong support for a cepheid like 
pulsation of this main mode, i.e. a fundamental 
radial mode. 



9.2. Theoretical periods of the the radial modes 
— Stellar models with "standard" mass loss rate. 



Lovy et al. (1984) and Schaller (1990) have calcu- 



lated the pulsation periods for the fundamental ra- 
dial mode and the first two overtones, from a linear 
non-adiabatic analysis. Schaller's analysis is based on 
the massive supergiant models from the grid of evo- 
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Table 6. Parameters of the highest amplitude variability modes according to equation (5) in various time intervals. 



Period [d] 


Frequency [c/d] 


HJD interval 


V Amp. A v 


[B - V] Amp. A B -v 


Av/Ab-v 


— (f>v — 4>B-V 


107.4 


0.00931 


2 447 850-2 448 400 


0.034 ±0.001 


0.017±0.001 


2.01 


-0.003 ±0.014 


156.4 


0.00639 


2 448 400-2 449 400 


0.059 ±0.003 


0.044±0.003 


1.34 


-0.014 ±0.018 


186.7 


0.00536 


2 447 850-2 450 100 


0.022 ±0.002 


0.014±0.002 


1.52 


0.071 ±0.038 




8000 
HJD-2'440'000 



Fig. 8. Three-frequency fit (0.00942 c/d, 0.00488 c/d and 
0.01232 c/d) for V and [B - V] {a and b respectively). 
Data are resampled in 3 days wide bin to homogenize the 
sampling. 



l ution ary tracks of massive stars by Schaller et al. 
( |l992| ). In these models, the dependence of the mass 
loss rate on L and T e g is that given by de Jager et 
al. ( |l988|) , called hereafter the "standard" mass loss 



rate. Schaller (g_990j) derived a theoretical PLC rela- 
tion, logP = -0.38M bo i - 3 log T cff ± 10.60, for the 
fundamental radial mode which predicts a period Pq 
= 140 d for V810 Cen. This value is valid for models 
cooler than logT e g = 4.1 and in the phase of helium 
burning. 

Stellar models with "high" mass loss rate. 

Meynet et al. ( 1994| ) have calculated a new grid of 
evolutionary stellar models, by adopting a "high" mass 
loss rate, i.e. a rate enhanced by a factor of two with 
respect to the "standard" mass loss rate. Various com- 
parisons with observations of high-mass stars support 
the high-mass loss grid of stellar models (see Maeder 



& Meynet, 1994). The linear non-adiabatic pulsa- 
tion code of Schaller ( |1992| ) has been applied to var- 
ious 25 M models, at the edge of the first redward 
evolution, close to V810 Cen position. The result of 



these calculations is presented in Table 7 for the lower 
(redward evolution) and upper (blueward evolution) 
tracks. A fundamental mode period of 157 d is ob- 
tained for the blueward evolution at logT c ff = 3.750 
and logL/L = 5.342, that is 0.004 dex away from 
values derived in section 4 for V810 Cen (logoff = 
3.78, logi/I/0 = 5.3). Due to a strong mass loss at 
the beginning of the blueward evolution, the mass has 
decreased to 20 M Q . 

9.3. Interpretation of the main modes 

— The observed main period at ^156 d. This ob- 
served period is in good agreement with the theo- 
retical prediction for models with "high" mass loss 
rate (157.2 d). Thus, the main observed period of 
V810 Cen (^156 d) can be explained by a pul- 
sation in the radial fundamental mode of a su- 
pergiant with an initial mass close to ~ 25 M Q . 
The evolutionary mass of V810 Cen must be ~ 20 M©. 
The observed second period at ^107 d. The pre- 
dicted periods of the first two radial overtones (Pi 
= 82.0 d, P 2 = 56.3 d, see Table 7) are much too 
small to explain this period. The period ratio of the 
observed two dominant modes is 0.69, while the pre- 
dicted ratios are P 1 /P = 0.52, P 2 /P Q = 0.36, P 2 /P 1 = 
0.69. This last value could suggest that V810 Cen pul- 
sates in the first and second radial overtones. However, 
in that case, the fundamental radial period would be 
^350 d. This period would require a very high lumi- 
nosity, Mboi — —9.5, which is difficult to postulate for 
a F8 la star. Since, in addition, this very large period 
is not observed, we conclude that the observed sec- 
ond main period cannot be a radial mode and, 
thus, a non-radial p-mode is the most natural 
explanation. 

The observed third period at ^187 d. This pe- 
riod is larger than the main mode at ~156 d. With 
the hypothesis that the main mode is the radial fun- 
damental pulsation, we have to conclude that the ob- 
served third period cannot either be a radial 
mode and, thus, a non-radial g-mode is to be 
postulated. Non-radial g-modes have already been 



sugges ted in supergiants by Maeder (1980) or de Jager 
et al. ( 1991 ), and are probably present in the ca se of 
yellow hypergiants like p Cas (Lobel et al., 1994 ). 



As noted above, the period ratio of the observed two 
dominant modes is 0.69, thus very close to the values of 
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the double-mode cepheids pulsating in the radial funda- 
mental mode and first overtone with Pi/Po m the 0.696- 
0.711 range (Balona, |l985| ). However, V810 Cen has a high 
mass, and therefore its Pi/Pq value is lower than the ratio 
observed for "low mass" double-modes cepheids (see table 
0). The similarity between the observed value of 0.69 and 
the 0.696-0.711 range for double-modes cepheids is at the 
origin of a wrong conclusion, i.e. V810 Cen is a double- 
mode cepheid in addition to its supergiant characteristics 
(Burki, |1994[ ). 

An alternative interpretation associating the ~187 d 
period to the radial fundamental has also been consid- 
ered; then the other modes would be non-radial p-modes. 
However, this alternative is unlikely because, if the funda- 
mental mode is excited, it should have the highest ampli- 
tude and this is not the case for the ^187 d period (the 
highest amplitude mode is the ~156 d period). Further- 
more, the (marginally significant) positive phase shift does 
not support that hypothesis neither and there are strong 
evidences for the radial nature of the ^156 d period. 

Our present conclusion, based on a very long-term pho- 
tometric monitoring, is that V810 Cen could be a super- 
giant star exhibiting the three types of pulsation modes : 
radial fundamental mode (main mode at ~156 d), 
non-radial p-mode (second mode at ^107 d) and non- 
radial g-mode (third mode at ~187 d). If this conclusion 
is correct, V810 Cen would be the first known case of a 
star showing such a diversity of pulsational characteristics. 

10. Conclusion 

Radial velocities, spectrophotometry and Geneva photom- 
etry data have been used to discuss the evolutionary sta- 
tus and pulsation modes of the yellow variable supergiant 
V810 Cen and its B-type companion. 

Radial velocities have been collected from the 
Coravel database as well as from the literature. They 
do not support the Stock 14 membership hypothesis of 
V810 Cen. To be consistent with the spectrophotometric 
data, the star has to be beyond Stock 14, at 3.3 to 3.5 kpc 
from the sun. 

IUE and visible spectrophotometry together with Ku- 
rucz atmosphere models have been used to estimate the 
physical parameters of the red supergiant : T c ft = 5 970 ± 
100 K, M v = -8.4,M bo i = -8.5, R = 420 R Q , spectral 
type F8 la. This star is about 3.3 mag brighter than its 
B0 III companion. Thus, the observed variability is due to 
the red supergiant. 

The derived temperature and luminosity place the star 
at the blue edge of the classical cepheids instability strip. 
The initial mass was ~ 25 Mq while the evolutionary mass 
ought to be close to ~ 20 M . 

The high-precision long-term photometric monitoring 
in the Geneva system reveals various types of variability : 

— A long-term increase of the luminosity from February 
1976 to July 1986. 



— A multimode pulsation behavior, with periods and am- 
plitudes variable with time (in the range 0.02 to 0.06 
mag). The main modes are : 

— A dominant ~156 d period, identified as the fun- 
damental radial mode. Its period value is in good 
agreement with the fundamental radial period de- 
rived from theoretical non-adiabatic linear pulsa- 
tion analysis. 

— The second mode with a period of ~107 d, which 
is most probably a non-radial mode, may be a p- 
mode. 

— The third mode with a period of ~185 d, which 
could be a non-radial g-mode. 

— Various other secondary modes have been detected 
during a time interval of high density measurements, 
at 89, 129, 167, 185 and 234 d. 

Are radial and non-radial p- and g-modes simultane- 
ously present in V810 Cen ? This is a quite fundamental 
question because this supergiant would then be the first 
known star exhibiting such a pulsational behavior. In our 
opinion, to answer this question, new hard and long-term 
work has to be done in the four following directions : i) 
New linear and non-linear non-adiabatic stability calcula- 
tions have to be done, based on well-adapted supergiant 
models; ii) High resolution (R~40 000) long-term spectro- 
scopic monitoring must be organized in order to follow 
the line profile variations; in) A long-term radial veloc- 
ity monitoring would help to determine the parameters 
of the orbit and of the components; iv) The photometric 
long-term monitoring must be continued in parallel to the 
spectroscopic survey. 
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Table 7. Periods and growth rates (rj) from linear non-adiabatic stability analysis of stellar model with "high" mass 
loss rate. Positive rj indicate unstable modes. 





log(T cff ) 




M 

M H 


[d] 


no 
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